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Neutrino diagnostics of ultra-high energy cosmic ray protons 
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The energy at which cosmic rays from extra-galactic sources begin to dominate over those from 
galactic sources is an important open question in astroparticle physics. A natural candidate is the 
energy at the 'ankle' in the approximately power-law energy spectrum which is indicative of a cross- 
over from a falling galactic component to a flatter extra-galactic component. The transition can 
occur without such flattening but this requires some degree of conspiracy of the spectral shapes and 
normalizations of the two components. Nevertheless it has been argued that extra-galactic sources 
of cosmic ray protons which undergo interactions on the CMB can reproduce the energy spectrum 
below the ankle if the cross-over energy is as low as the 'second knee' in the spectrum. This low 
cross-over model is constrained by direct measurements by the Pierre Auger Observatory which 
indicate a heavier composition at these energies. We demonstrate that upper limits on the cosmic 
diffuse neutrino flux provide a complementary constraint on the proton fraction in ultra-high energy 
extra-galactic cosmic rays and forthcoming data from IceCube will provide a definitive test of this 
model. 



PACS numbers: 98.70.Sa, 95.85.Ry 



I. INTRODUCTION 



The flux of cosmic rays (CRs) falls as an approxi- 
mate power-law in energy, dN/dE cx E^ 1 , with 7 ~ 2.8 
from about 1 GeV up to the 'knee' in the spectrum at 
~ 3 x 10 6 GeV where it steepens to 7 ~ 3; it then steep- 
ens further to 7 ~ 3.2 at ~ 5 x 10 8 GeV (the 'second 
knee') and flattens back to 7 ~ 2.8 at the 'ankle' at 
~ 3 x 10 9 GeV [HQ. A long-standing open question is 
the transition point between dominance by galactic and 
extra-galactic sources in the spectrum of ultra-high en- 
ergy (UHE) CRs. The transition between the two com- 
ponents ought to be accompanied by the appearance of 
spectral features, e.g. two power-law contributions would 
naturally produce a flattening in the spectrum if the 
harder component dominates at lower energies. Hence, 
the ankle seems to be a natural candidate for this tran- 
sition [H, 0, [f| (for discussions see Refs. [1, 0|)- 

It has been argued however that the cross-over of the 
two contributions might also happen at lower energies 
without a corresponding flattening of the spectrum or 
even at a spectral steepening like the second knee 0. 
Such a transition would appear to require considerable 
fine-tuning of the shape and normalization of the two 
spectra. It has been argued however that the spectral 
shape required for the extra-galactic component can de- 
velop naturally during the propagation of extra-galactic 
protons in the cosmic microwave background (CMB) over 
cosmological distances. At an energy of ~ 3 x 10 8 GeV 
the energy losses due to e + e~ pair production and cosmic 
expansion are roughly equal and thus produce a steep- 
ening of an initially featureless power-law injection spec- 
trum. In this model the ankle is formed as a result of 
the dip due to e + e~ pair production together with a pile- 
up of protons below the Greisen-Zatsepin-Kuzmin (GZK) 
suppression which sets in at about 6 x 10 10 GeV due to 



pion photoproduction [^, [Io| . 

An obvious test of this 'dip-transition' model Q is 
whether the flux is indeed proton dominated at en- 
ergies above 10 9 GeV [ll|. Recent measurements of 
the elongation rate of UHE CR air showers by the 
Pierre Auger Observatory [12j indicate however an 
increasingly heavier composition in the energy range 
- 2 x 10 9 - 4 x 10 10 GeV [H[, in conflict with this ex- 
pectation. More elaborate models [3, EH EH m which 
the the extra-galactic CRs have a mixed composition and 
the cross-over occurs between the second knee and the 
ankle are consistent with the data [Tt], Ell , but lack the 
simplicity of the dip-transition model. However, this il- 
lustrates that even precision measurements of the chemi- 
cal composition and spectrum of CRs may not be enough 
to resolve conclusively just where the galactic-to-extra- 
galactic transition occurs. 

An important clue in resolving this puzzle would be the 
identification of the CR accelerators themselves. Indeed, 
the Auger data show a correlation (at > 99% C.L.) be- 
tween the arrival directions of CRs with energy above 
~ 6 x 10 10 GeV and active galactic nuclei within 
~ 75 Mpc pjj [2(| • However, this has only deepened the 
puzzle, indicating the possible return of a lighter compo- 
sition beyond 4 x 10 10 GeV, since the primaries appear 
not to have been deflected significantly by intervening 
inter-galactic and galactic magnetic fields and therefore 
are more likely to be protons rather than heavy nuclei. 
The cosmic ray sources ought also to be emitting other 
messenger particles that are produced during the accel- 
eration process, in particular gamma rays and neutri- 
nos, and these might provide further constraints. In- 
deed, it has been observed [2l[ that certain low cross-over 
models which require very powerful sources are already 
marginally excluded by experimental upper bounds on 
diffuse neutrino fluxes. However, a general exclusion of 
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such models based on neutrino bounds alone has not yet 
been possible, given the uncertainties inherent in the ac- 
celeration mechanism as well as the possible cosmic evo- 
lution of the sources. 

We update and improve this observation [2l| in sev- 
eral respects. The bound on diffuse neutrino fluxes be- 
low 10 9 GeV has improved by a factor of 4 between the 
AMANDA-BIO [13 and the AMANDA-II [H H mea- 
surements and will improve further by an order of mag- 
nitude after 1 year of observation with IceCube which is 
currently under construction at the south Pole [25[ . In 
addition, the hybrid measurements of the Pierre Auger 
Observatory have reduced the uncertainties in the abso- 
lute normalization of the UHE CR flux [2(| and these 
measurements disagree with both the AGASA data [27| 
and the HiRes data 0, [29[ , which was used in the ear- 
lier analysis [2l[ (see also Ref. [30]). Furthermore, we 
will present a method to derive differential upper limits 
on the proton fraction of UHE CRs from experimental 
bounds on the diffuse neutrino flux. 

The outline of this paper is as follows. We start in 
Sect. IH by der iving the luminosity relation between neu- 
trons and neutrinos in optically and magnetically thin 
sources. In Sect. [Hi] we discuss the Boltzmann equations 
that govern the propagation of UHE CRs and neutrinos. 
These results are used in Sect. IIVI to derive upper limits 
on the extra-galactic proton contribution to UHE CRs 
from candidate sources such as blazars. We summarise 
our findings in Sect. [V] and provide derivations of the 
results referred to in the main text in the Appendices. 



II. EXTRA-GALACTIC NEUTRINO 
PRODUCTION 



Cosmic rays originating in cosmic accelerators will typ- 
ically be accompanied by gamma-rays and neutrinos. 
This is a consequence of the inelastic hadronic processes 
which are involved in their production mechanism [3l| . 
Candidate sources for the highest energy CRs such as 
blazers or gamma ray bursts are expected to accelerate 
charged particles by the lst-order Fermi process or 'diffu- 
sive shock acceleration' in which they repeatedly scatter 
off a propagating plasma shock front (see Ref. [32|]). (It is 
also possible that UHE CRs are accelerated by the 2nd- 
order Fermi process, e.g. in the extended lobes of radio 
galaxies [HI].) For efficient acceleration up to the highest 
observed energies this process has to be repeated many 
times hence it is essential that a magnetic field confines 
the charged particles for sufficiently long in the vicinity 
of the shock front. 

Accelerated electrons and other light charged parti- 
cles will lose their energy in the magnetic field due to 
synchroton radiation. The resulting photons provide a 
target for protons and heavier nuclei to undergo meson 
photo-production and photo-disintegration, respectively. 



A neutron produced in this process may escape out of 
the magnetically confined source before it /3-decays into 
a CR proton. To maintain the efficiency of the accelera- 
tion mechanism, such interactions need to be much less 
frequent than the acceleration cycle. Charged and neu- 
tral pions decay further into high energy neutrinos and 
gamma rays which are not confined by the magnetic field 
and can be emitted from the source. 

Here we focus on cosmic proton accelerators. De- 
pending on the relative ambient gas and photon den- 
sities, charged pion production may proceed either 
through inelastic pp scattering |34j |. or photopion pro- 
duction predominantly through the resonant process 
P7 — > A + — > nir + [35| . The subsequent decay chain 
7r + — » — * e + v e V ti v il produces high energy neutrinos 
with an initial flavor ratio v e : : v T of 1 : 2 : 0, which is ex- 
pected to evolve through oscillations into approxim ately 
1:1:1 after propagation over cosmological distances [36| . 

The relative luminosity of neutrons and neutrinos de- 
pends on various details of the source [3l[. Since the 
interaction time scale of protons and neutrons is simi- 
lar, the relative luminosity of neutrinos will tend to be 
higher in sources which are optically thick to (p/n)j in- 
teractions. The life-time, geometry, bulk motion, and 
magnetic field of the source limits the maximal energy 
of the emitted protons [37(. Also important are cool- 
ing processes of the muons and pions in the source prior 
to decay, such as adiabatic losses due to the expansion 
of the source, inverse Compton scattering, Bethe-Heitler 
pair production, and synchrotron radiation losses [3lT ] . 

Our working hypothesis in the following can be ex- 
pressed in terms of the time scales of pry and pp inter- 
actions (t pp and r P7 ), the decay lifetimes of neutron, 
pion and muon (r„, tv and r^), the cooling time scale of 
charged particles (r coo i), the characteristic cycle time of 
confinement (r cyc i c ), and the total confinement time scale 



(^conf)- (i) 7p7 Tpp: (^v 7"cycle ^ ^~ni (Hi) 7"cycle 

(iv) t v1 < T conf ; (v) T n/7r /^ < 7"cool- 



< 



'P7, 



Conditions (i) to (Hi) ensure that photopion produc- 
tion of neutrons and subsequent diffusion out of the mag- 
netic confinement region is the dominant process of cos- 
mic proton production in an optically thin source. Con- 
ditions (Hi) and (iv) permit efficient cosmic ray accel- 
eration with a smooth spectrum across a wide range of 
energy. Finally, condition (v) ensures that cooling pro- 
cesses of secondary neutrons, pions, and muons have neg- 
ligible influence on the energy distribution of the emitted 
neutrons and neutrinos. 

Within this working hypothesis the relative number 
and energy of the neutrinos and neutrons depend only 
on the kinematics of photo-hadronic interaction, which 
implies approximate equipartition of the decaying pion's 
energy between the neutrinos and the electron, and the 
relative radiation density in the source. On average, each 
proton-photon interaction will produce 77 neutrinos per 
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FIG. 1: The values of e = (E„) / (E n ) and T) = <AU/(iV„) as 
calculated with the SOPHIA Monte Carlo code [|cj] for the 
blazar photon target spectrum shown in Fig. 1 of Ref. [jjjj. 
For comparison we show the values at the A + resonance used 
previously in Ref. [2l| . 



neutron with relative energy e per neutrino, i.e. 



V 



(Nu) 

(Nn) 



and 



(En) 



(1) 



In the following, we will consider a (hypothetical) source, 
where photo-pion interaction proceed exclusively via the 
A + resonance with fixed values rj = 3 and e = 0.07 [2l| . 
and the blazar flaring state model [38| that has been pre- 
viously discussed in Ref. [39j . 

For a given photon target spectrum in the source the 
values of n and e can be obtained using the SOPHIA 
Monte Carlo code [4(| for photo-nucleon interactions and 
are shown in Fig. [T] for the adopted blazar spectrum. For 
a fixed proton energy we sample the scattering angle and 
photon energy in the center-of-mass frame, weighted by 
the total P7 cross section and the photon distribution, 
i.e. the integrand of Eq. (|Aip . These parameters are 
then input to SOPHIA to generate a pj interaction and 
the process is repeated 10 5 times to derive the average 
values for the neutrino/neutron number and energy. 



The neutrino emissivity of flavor i is then given by: 



AE V 



■C Vi {z,E v 



AE n 



(2) 



Assuming flavor universality as well as e ~ E v /E n ~ 
AE V / AE n and rj ~ N a \\ v /N n we arrive at the neutrino 



source luminosity (per co- moving volume) : 



Caa.v{z,E v ) ~ - C n (z,E v /e) 
e 



(3) 



Note, that the relation §5§ derived for optically thin 
sources can be regarded as a lower limit on the neutrino 
luminosity as long as energy-loss processes in the source 
are negligible. We will illustrate this below for the vari- 
ous conditions of our working hypothesis one by one. 

Pion production may also proceed predominantly via 
inelastic pp scattering if r p7 > t pp (cf. (i)), in the source. 
For the process pp — » NN + pions, hadronic event gen- 
erators indicate an inelasticity ~ 0.6 where N is a final 
state nucleon [4l|, H2]- Assuming that 2/3 of the final 
state pions are charged we estimate that the average en- 
ergy deposit into neutrinos is about rje ~ 1/4 with r\ 3> 3 
due to the large multiplicity of secondary pions in inelas- 
tic pp collisions. For a neutron luminosity C n oc E^ 1 
with 7^2 (typical for relativistic plasma shocks) this 
will somewhat increase the neutrino flux relative to the 
neutrons (cf. Eq. ©). 

If neutrons decay within the source, i.e. T cyc \ c > r n 
(cf. (ii)), and the sources are optically thick, i.e. T cyc \ c 
t pi ( c f- then the neutrino to neutron ratio would 

also be enhanced due to neutron re-conversion into pro- 
tons through /3-decay and nj interactions, as well as 
simultaneous production of additional neutrinos. The 
optical thickness of the source depends on the intensity 
of the radiation field as well as the characteristic size 
R of the acceleration region. In particular, gamma-ray 
bursts, possible sources of UHE CRs, are optically thick 
for proton energies larger than 10 7 GeV given standard 
parameters [43l[44| so we do not consider direct neutrino 
emission from such sources in the following (conservative) 
analysis. 

However, if the target photon field gets too thin so 
that T C ycie -C t P7 (cf. (Hi)), then the source would rather 
accelerate protons to higher energy than emit neutrons 
through photopion interactions. Whereas this does not 
change the ratio (J3J) and our analysis remains conserva- 
tive, both neutron and neutrino emission become ineffi- 
cient. This is important for neutrino energies below a few 
10 8 GeV in blazars [39[, hence we will not use experimen- 
tal neutrino bounds below 10 8 GeV for these candidate 
sources. However, one can envision cosmic ray sources 
in which there is substantial leakage of protons from the 
vicinity of the shock wave, which however remain trapped 
in the source magnetic field until they interact with the 
photon background, producing neutrons and neutrinos. 
Thus, in our general analysis of optically thin sources we 
extend the energy range down to 10 6 GeV. 

Some care has to be taken if cooling processes of sec- 
ondary pions and muons in the source environment are 
in fact important (cf. (v)) [U Hf|. Synchrotron radi- 
ation, inverse Compton (IC) scattering, Bethc- Heftier 
(BH) pair production, or adiabatic losses can reduce the 
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energy of pions and muons before decay and hence de- 
crease the neutrino luminosity at high energies. Inverse 
Compton cooling is the most relevant process for pho- 
tointeractions of pions and muons [3l| and its time scale 
can be related to the synchrotron loss time via the rela- 
tive co-moving energy densities of photons and the mag- 



netic field, 



4C _ 



sync sync 



T* ync /£- Here, we have 



introduced the equipartition parameter £ which can ex- 
ceed 1 in baryon-loaded flows due to the low radiation 
efficiency of relativistic protons 31]. However, Klein- 
Nishina corrections at high pion and muon energies re- 
duce the efficiency of IC cooling, which might effectively 
lead to £ < 1 [4q . Hence, we will estimate the time-scale 
of secondary particle cooling by r* ync . 

Synchrotron losses of secondary pions and muons in the 
background magnetic field are important if their lifetime 
7*T^ ec is larger than the energy loss time T* ync , given by 1 



sync 



rnZ 



16tt a 2 E+B 2 



(4) 



With muon and pion lifetimes of t^ cc = 2.2 /is and 
26 ns, respectively, this translates into a critial en- 
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above which our relation between neutron and neutrino 
luminosities is unreliable. 

We use neutrino bounds from AMANDA II [H, [I| 
and Auger as well as the projected sensitivity of Ice- 
Cube [43], which extend up to energies of a few 10 10 GeV 
(in the observatory frame) and, after red-shifting, corre- 
spond to neutrino emission in the source up to energies of 
about 10 11 GeV and pion/muon production up to a few 
10 11 GeV (since the neutrinos carry only about 1/4 of 
the pion/muon energy). Equation (O shows that cooling 
processes in blazar jet environments with typical boost 
factor T ~ 10 15 and B < 10 G [HIH may become im- 
portant for neutrinos from muon decay only at energies 
of about 10 10 GeV. However, it has been shown [39[ that 
typically blazar sources become optically thick at these 
energies which leads to an increase of the neutrino to pro- 
ton ratio. In contrast, GRBs typically have a much larger 
magnetic field of order B ~ 10 12 G and a boost factor of 
r ~ 10 25 . Hence, cooling of secondaries in GRBs can- 
not be neglected, but as we have noted already they are 



not optically thin sources at UHE and we do not there- 
fore consider (direct) neutrino emission from them in the 
following (conservative) discussion. 

The ratio of the observable flux of photo-hadronic neu- 
trinos and CRs today can be calculated from the relation 
([3]) using a set of Boltzmann equations that take into 
account the interactions with the cosmic photon back- 
ground and the dilution and red-shift via cosmic expan- 
sion. We will discuss these equations and their solutions 
in the next section. 



III. PROPAGATION OF COSMIC RAYS AND 
NEUTRINOS 



For a spatially homogeneous distribution of cosmic 
sources, emitting UHE particles of type i, the co-moving 
number density Yj is governed by a set of 1-dimensional 
Boltzmann equations of the form 

Yi = d E {HEYi) + d E (biYi) 

I dE j 1]l Y ] + d , (6) 
3 J 

together with the Friedman-Lemaitre equations describ- 
ing the cosmic expansion rate H{z) as a function of the 
redshift z. Non-relativistic and non- interacting matter 
obeys the equation Y = 0. The first and second term 
in the r.h.s. of Eq. © describe continuous energy losses 
(CEL) due to red-shift and e + — e~ pair production on 
the cosmic photon backgrounds, respectively. The third 
and fourth terms describe more general interactions in- 
volving particle losses (i — > anything) with interaction 
rate Ti, and particle generation of the form j — > i. 2 The 
last term on the r.h.s. of Eq. (J6j) , £,-, corresponds to the 
luminosity density per co-moving volume of sources emit- 
ting CRs of type i. In Appendices E] [B] and [C] we pro- 
vide more details regarding the quantities appearing in 
Eq. ©. 

Note, that the Boltzmann equations (J6j) do not take 
into account the deflection of charged CRs during their 
propagation through inter-galactic and galactic magnetic 
fields. In fact, if synchrotron radiation during propaga- 
tion is negligible and the source distribution is homoge- 
nous, Eq. ((6]) provides a good approximation of the spec- 
tral evolution even for CRs having small rigidity which 
suffer large deflections [50| . However, magnetic inhomo- 
geneities on small scales will suppress the spectrum of 
CRs with Larmor radius l\, < id where id is the char- 
acteristic distance between sources. It has been shown 
that for typical inter-galactic magnetic fields of strength 



We work throughout in natural Heaviside-Lorentz units with h = 

c = eo = MO = 1) a = e 2 /(47r) ~ 1/137 and 1 G ~ 1.95 X 2 We will illustrate in Appendix 1X1 that CEL processes are an ap- 
10 _2 eV 2 . proximate formulation of the process i — » i. 
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FIG. 2: Integrated and differential upper limits on the proton contribution in UHE cosmic rays. The thin lines show a 
sample of proton and neutrino test spectra (see Eqs. (|3]) and $12$ ) used in the analysis. Upper panels: Limits derived from 
the AMANDA-II bound on diffuse neutrinos [§J, E 2 J t3Xv < 2.7 x 10" 7 GeV cm" 2 s _1 sr" 1 (90% C.L.) from 2 x 10 5 GeV to 
10 9 GeV, assuming that the sources are optically thin with 77 = 3 and e = 0.07 (c/. Ref. [IJ]). Lower panels: Limits derived 
from the Auger bound on UHE tau neutrinos [3], E 2 J^ +i , T < 1.3 x 10" 7 GeV cm" 2 s" 1 sr" 1 (90% C.L.) from 2 x 10 8 GeV to 



2 x 10 GeV, assuming that the sources are blazars. We adopt a neutrino flavor ratio v e 



:v T of 1 : 1 : 1 at Earth to exploit 



this bound for other flavours and to extract the differential limit we use the (relative) value of the exposure quoted in Ref. |59|. 



~ 1 nG and coherence length of ~ 1 Mpc, the diffu- 
sive propagation of CR protons will start to affect the 
spectrum below about 10 9 GeV if £ d ~ 50 Mpc De- 
pending on the diffusion regime, this can suppress the 
proton flux at 10 8 GeV by a factor of 3 to 100. We will 
return to this point in the discussion of our results. 

We adopt the usual concordance cosmology of 



a fiat universe dominated by a cosmological constant 
with £1a ~ 0.7, the rest being cold dark matter 
with f2 m ~ 0.3. The Hubble parameter is given by 
H 2 {z) = Hq (£! m (l + z) 3 + £1a), normalised to its value 
today of 70 kms" 1 Mpc" 1 . The time-dependence of the 
red-shift can be expressed via dz = — dt (1 + z)H . 

The cosmological evolution of the source density per 
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FIG. 3: As in Fig. [2] but now showing the estimated sensitivity to the proton component in UHE CRs from the IceCube 
reach after one year of observation. The IceCube sensitivity E 2 J^\\ v ~ 2.4 x 10 -8 GeV cm -2 s _1 sr _1 on UHE neutrinos above 
10 8 GeV is estimated by an extrapolation to 10 10 GeV motivated by the analysis of Ref. [4?| . Upper panels: Sensitivity 
to the proton component exploiting both cosmogenic neutrinos and neutrinos from optically thin sources. Lower panels: 
Sensitivity exploiting cosmogenic neutrinos alone in models where the CR sources evolve strongly with redshift; we show the 
case of strong cosmological evolution of the proton sources according to Hgrb ([9] left panel) and Wagn (1101 right panel). 



co-moving volume is parameterized as: 

jC i (z,E)=n(z)C i (0,E), (7) 

where the source luminosity per co-moving volume is as- 
sumed to follow the star formation rate (SFR). (Note 
that the dilution of the source density due to the Hubble 
expansion is taken care of since C is the comoving den- 
sity, hence for no evolution we would simply have Ti = 1.) 



Following the recent compilation [5^, HH we adopt 

f(l + z) 3 - 4 z<l, 
HsfrO) - < iVi(l + z)- - 3 Kz<4, (8) 
[iViJV 4 (l + z)- 3 - 5 z>4, 

with appropriate normalization factors, N% = 2 3 7 and 
-/V 4 = 5 3 2 (see the right panel of Fig. [5]). 
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Some candidate sources of UHE CRs may have a 
stronger evolution than the star formation rate ||5J); 
this is particularly important for cosmogenic neutri- 
nos [H, HH, [56[ . For gamma-ray bursts we adopt [HiJ 



Hgkb(z) = (1 



l n SFR (z) 



(9) 



Active galactic nuclei may have a similarly strong evolu- 
tion. Following Refs.[H, H3| we take 



(1 + z) 5 - z<1.7, 
W\,;n(--1 = {N U7 1.7<z<2.7, 
JVi.7 N^ 2 7 7 ~ Z) z > 2.7 , 



(10) 



with iVi.7 = 2.7 5 and iV 2 . 7 = 10 43 . We will discuss the 
impact of these strong evolution scenarios on cosmogenic 
neutrinos and the proton fraction at very high energies. 

In Appendix [D] we derive a general solution to the dif- 
ferential equation © . The neutrino flux from the sources 
is obtained simply by integrating Eqs. (|D1[) and (|D2|) : 



ME) = dzj±-r£ u (z, (1 + z)E) . (11) 



To obtain the cosmogenic neutrino flux produced during 
propagation through cosmic radiation backgrounds we 
need to solve the full Boltzmann equations numerically. 



with an exponential energy cut-off i? max that we vary 
between 10 8 GeV and 10 12 GeV with a logarithmic step- 
size of log 10 E = 0.25. Unless otherwise stated, we use 
Eqs. ([7]) and © for the source evolution. Each neutron 
test luminosity (|12p is related to a neutrino luminosity 
by the ratio |3]). After propagation using Eqs. © we 
normalize the accumulated contribution of extra-galactic 
and cosmogenic neutrinos to the limit on the diffuse neu- 
trino flux from AMANDA II [2! and to the limit on UHE 
tau neutrinos from Auger [4q|. This can be done in two 
possible ways. 

For an integrated upper limit (or sensitivity) we max- 
imise the individual neutrino flux normalization to sat- 
urate the integrated experimental neutrino flux bounds. 
This approach is shown in the left panels of Fig. [21 The 
envelope of the corresponding proton flux gives the inte- 
grated upper limit on the proton fraction. However, the 
application of this limit requires that the flux of cosmic 
neutrinos is close to an E~ 2 spectrum in the quoted in- 
terval [E-,E+\. Note, that for CR protons above a few 
times 10 s GeV this upper limit becomes trivial, i.e. it lies 
above the observed CR flux. 



A differential upper limit can be obtained using the 
effective area A e g of AMANDA-II, provided in Ref. 
and the acceptance (oc A c g) of Auger quoted in Ref. 
The normalisation of the neutrino test fluxes J* est is fixed 
by the integrated limit J" 1 * and the equation 



24], 
59 j. 



IV. BOUNDS ON THE EXTRA-GALACTIC 
PROTON FRACTION 



dEA eS J^ 



&EA eS Jl cst 



(13) 



The luminosity relation between neutrinos and 
CRs in optically thin sources was first used [35j to de- 
rive an upper limit on the diffuse neutrino background 
from UHE CR observations. This argument was modi- 
fied and extended in further studies which examined if 
the upper limit on the neutrino flux could in principle 
be exceeded [U, [58| . With the advent of km 3 -scale neu- 
trino telescopes like IceCube these limits will soon be 
tested dj. 

As we have outlined in Section |Hj the relative neu- 
trino luminosity ([3]) of optically thin sources with photo- 
hadronic neutrino production (conditions (i) to (Hi)) can 
be regarded as a lower limit on the neutrino flux, if cool- 
ing processes of secondary pions and muons are negligible 
(iv) and proton diffusion out of the confinement region is 
suppressed (v). We translate this conservative expecta- 
tion into an upper limit on the extra-galactic proton frac- 
tion in UHE CRs, exploiting experimental upper bounds 
on the diffuse high-energy neutrino flux. 

For this procedure it is convenient to introduce test 
functions of the neutron source luminosity of the form 



C t : st (0,E) = C o 



E 



E n 



-1 
exp 



E 



EL 



(12) 



with the appropriate integration limits E± according to 
the experimental integrated neutrino flux bound. This 
second approach is shown in the right panels of Fig. [5] 
Again, the envelope of the corresponding proton fluxes 
provides the differential upper limit on the proton frac- 
tion. Note, that this differential bound is more general 
and does not make the usual assumption that the neu- 
trino energy spectrum is oc E~ 2 , but is in general weaker 
than the integrated limit by up to a factor of 10. 

One year of observation with the fully deployed Ice- 
Cube detector should increase the sensitivity to diffuse 
neutrino fluxes by about an order of magnitude below 
10 8 GeV. It is conceivable that this sensitivity can be ex- 
tended to UHE neutrino energies of about 10 10 GeV [i^ . 
The corresponding (integrated) sensitivities on the pro- 
ton fraction in UHE CRs are shown in the upper panels 
of Fig. [3] We see that 1 year of observation is sufficient to 
place stringent bounds on the extra-galactic proton frac- 
tion in CRs below 10 10 GeV and thus provide a definitive 
test of the 'dip-transition' model Q • 

As noted earlier, the interplay of a random walk of 
cosmic protons in magnetic fields together with an in- 
homogenous distribution of their sources may result in 
suppression of the observed proton flux below about 
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FIG. 4: Summary of differential upper limits (thick lines 
1&2) and integrated sensitivities (thin lines 3-6) on the 
extra-galactic proton contribution compared to CR data from 
HiRes-I/II (recalibrated) and the Pierre Auger Observatory 
(see Figs. [5] and |3] for details). The limits and sensitivities 
are derived from the A + approximation (1&3), from blazar 
emission (2&4) and from cosmogenic neutrinos with a strong 
source evolution Hagn (5) and Hgrb (6), respectively. 



CONCLUSIONS 



Exactly where the transition occurs between the galac- 
tic and extra-galactic components of UHE cosmic rays is 
presently an open question and identifying this energy 
would provide important clues to unravelling their ori- 
gin. The acceleration of CRs in extra-galactic sources 
would be accompanied by the emission of high energy 
neutrinos, thus an upper limit on the proton flux can 
in principle be inferred from experimental bounds on the 
extra-galactic UHE neutrino flux. We have demonstrated 
this, focussing on blazar jets which we have argued are 
optically thin sources, so photo-hadron interactions gen- 
erate a certain minimum flux of UHE cosmic neutrinos 
as long as proton diffusion and secondary particle cool- 
ing are negligible. We have shown that this argument is 
conservative in that a higher neutrino flux would be ex- 
pected if any of our assumptions are relaxed, e.g. if the 
sources of UHE CRs are optically thick. 

Our main results are summarized in Fig. 3J The 
AMANDA-II bound on diffuse neutrinos already con- 
strains the extra-galactic proton contribution in CRs at 
energies below a few times 10 8 GeV and just 1 year of ob- 
servation with IceCube will provide the necessary sensi- 
tivity up to a few times 10 10 GeV. If the number density 
of extra-galactic CR sources evolves strongly with red- 
shift then the detection of cosmogenic neutrinos alone 
may enable bounds to be placed on the proton flux at 
trans-GZK energies. 



10 9 GeV [5l|. Hence, we expect that our upper limit 
on the proton flux in this energy region is rather con- 
servative. In addition, a source evolution much stronger 
than the star formation rate ([5]) would predict a much 
higher flux of neutrinos - cosmogenic as well as from the 
sources - whereas the spectrum of protons is fixed by the 
contribution of local sources within the 'GZK horizon'. 

The flux of cosmogenic neutrinos, which dominates the 
total neutrino flux above 10 10 GeV, can also in princi- 
ple be exploited to obtain bounds on the proton frac- 
tion independent of the nature of the sources [55| . How- 
ever, since the CR flux is steeply falling as ~ E~ 3 , to 
constrain the contribution of trans-GZK protons would 
require even more sensitive experiments than IceCube 
and/or a significantly stronger source evolution with red- 
shift than we believe is plausible. In the lower panels of 
Fig. [3] we estimate the sensitivity of IceCube to constrain 
the proton fraction using only cosmogenic neutrinos for 
two strong evolution scenarios (see Eqs. © and (110p ). 
Again, the IceCube sensitivity shown to UHE neutrinos 
up to 10 10 GeV is extrapolated from lower energies, mo- 
tivated by the analysis of Ref. (47[. We observe that 1 
year of observation at IceCube would be sufficient to con- 
strain trans-GZK protons in UHE CRs only if the source 
evolution is sufficiently strong: oc (1 + z) 5 or steeper. 
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APPENDIX A: PHOTON BACKGROUND 
INTERACTIONS 



The propagation of UHE cosmic nuclei is affected 
by photo-hadronic interactions on cosmic photon back- 
grounds. For UHE protons the dominant interactions 
that determine the spectrum occur on the CMB, but 
the cosmic infra-red background (CIB) is also impor- 
tant for the generation of cosmogenic neutrinos. For the 
spectrum of the latter we use the recent compilation of 
Ref. 
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E [GcV] red-shift z 

FIG. 5: Left Panel: The interaction and decay rates appea ring in the Boltzmann equations for the CMB and CIB [6(j at 
z — 0. Right Panel: Star formation rate (Eq. (|8]) from Ref. J53|] ) and our approximation of the CIB number density scaling 
with redshift (|C2|I , For comparison, we also show the scaling behaviour of the CMB number density oc (1 + zf . 



The angular-averaged (differential) interaction rate 1\ 
(jij) appearing on the r.h.s. of Eq. ^ is defined as 

l 

Ti(z, Ei) = X - J dc e J de (1 - 0c e )n^(z, , (Al) 
-l 

7ii (z, Ek, Ej) = Ti(z, Ek) ^(E^Ef , (A2) 

where n 7 (z, e) is the energy distribution of background 
photons at redshift z and dNij/dEj is the angular- 
averaged distribution of particles j after interaction. For 
photo-hadronic interactions this distribution can be de- 
termined using the Monte Carlo package SOPHIA [io| . 
The factor 1 — [3ce takes into account the relative motion 
of photons and the nucleus, i.e. the Doppler shift of the 
photon density. 

We assume that the photon background has the adia- 
batic scaling behaviour: 

n 7 (z, e) = (1 + zf n 7 (0, e/(l + z)) . (A3) 

This is exact for the CMB (following from 
Yy = dsiHEY^) and Yy oc a 3 n 7 ), but not so for 
the CIB. However, the dominant opacity for proton 
propagation is provided by the CMB. The scaling 
behaviour Eq. (|A3|) translates into the following scaling 
of the quantities Ti and "fij , 

Ti(z, Ei) = (1 + zf T^O, (1 + z)Ei) , (A4) 

7y («, E^Ej) = (1 + zf 7 y (0, (1 + z)E h (1 + z)Ej) . 

(A5) 



For the adopted scaling of the CIB see Appendix [Cl 

If all interactions can be described as a CEL pro- 
cess the differential equation (|D2[) is considerably sim- 
plified and can be solved in a closed form as we will 
show later. In general, any transition i — > i which can 
be approximated as "/ U (E,E') « S(E — E' — AE)Ti(E) 
with AE/E <C 1 can be replaced in the Boltzmann equa- 
tions ([!]) as 

-T(E)Y t (E)+ J dE' 7u (E',E)Yi(E') 

-^d E (biY), (A6) 

with bi = AETi —E. The production of electron- 
positron pairs in the photon background with a small 
energy loss is usually approximated as a CEL process. 
Here we follow the standard approach of Ref. [61( to cal- 
culate the proton energy losses in the photon background. 
Again, the computation of the quantity b at various red- 
shift is significantly simplified if we assume an adiabati- 
cally scaling background photon density as for the CMB. 
The scaling behaviour of b and its derivative — d E b is 
then 

h(z, Ei) = (1 + zf 6,(0, (1 + z)Ei) , (A7) 
0i{z, Ei) = (1 + zf pi{0, (1 + z)Ei) . (A8) 

As before, for the scaling of the infra-red background see 
Appendix [C] In the left panel of Fig. ^ we show the 
quantities 6P air /E, /3 pan , T p and Hq for comparison. 
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APPENDIX B: NEUTRON DECAY 

For neutrons produced by proton interactions on the 
photon background, the Boltzmann equations have to in- 
clude a decay term. Since the proton and neutron are 
much heavier than the electron and electron anti-neutrino 
we can safely assume that the proton is at rest in the cen- 
ter of mass frame and Ep + E* « m n — m p = Am. The 
energy distribution of electron and electron anti-neutrino 
in the nucleon's rest frame can then be approximated as 



dN e 

dE* 
dN p 

dE* 



oc 



oc 



E* 2 -m 2 e E*(Am- E*f 



(Bl) 



E 



*\2 



m 2 e E* z (Am-E*), (B2) 



according to the phase space density of electron and 
anti-neutrino with dN ~ d 3 p e d 3 pp c and the constraint 
Am = E e + Ep e . The angular-averaged distribution in 
the lab frame is 



where £cib is the co-moving luminosity density of the 
sources and we neglect absorption effects other than ex- 
pansion. We assume that this follows the star formation 
rate: Ccib(z,E) oc Hsfr(z) £cib(0, E). We can then 
derive the bolometric evolution as 



Ncmjz) 
JVcib(O) 



= (1- 



,r 3 dz'?w(ff(^)(i+^)) 

f™dz>H S FR/(H(z>)(l + z>)) ' 



(C2) 



where Ncib(z) is the number of infrared-optical photons 
per proper volume at red-shift z. For comparison, the 
CMB evolves as Womb (^)/-Wcmb (0) = (1 + zf. To sim- 
plify the numerical evaluation we approximate the evo- 
lution with rcdshift as 



n cm (z,e) 



1 N cm (z) 

T— -r T tttt «cib(0, e/(l 

1 + z A^cibIO) 



-z)). (C3) 



dN 
dE 



dew 



dE* dN 
~dEdE*~ ' 



(B3) 



with E = 7„-B* — JnPnCeP*- The energy of the emerging 
proton is approximately j n m p = m p E n /m n ~ E n . 

Neutron decay contributes then in the Boltzmann 
equation in the quantities 

m 

dE 



and 



In 



= r " cc -r^- (B4) 



for i = D e ,e~. For neutrons with energy less than 
10 11 GeV the decay length is always smaller than the in- 
teraction length in the photon backgrounds (c/. left panel 
of Fig. [5]). In this case it is convenient to approximate 
the production of neutrons as 



r 



pe 



peff p7r 



= 

pe 



dEn^p n 

dE n T 



pn 



dN e - 
dE 
dNp e 
dE 



PP PP P n ' 



(B5) 
(B6) 
(B7) 



We have used this approximation in all our calculations. 



APPENDIX C: EVOLUTION OF THE 
INFRARED-OPTICAL BACKGROUND 



The CIB spectrum has been studied and tabulated in 
Ref. [60j for red-shifts up to z = 2. The red-shift depen- 
dence is given by 



n C m(z,(l + z)E) 

OC 

= (1 + z) 2 / dz' 



H(z') 



-CIB 



(z',(l + z')E), (CI) 



The redshift scaling of the quantities 7^, T^, bi and 
fa for the CIB is then obtained from the corre- 
sponding scaling given for the CMB in Eqs. (|A4^IA5j) 
and (|A7| / IA8|) . by multiplying the r.h.s. by a factor 
Ncib(z)/Ncib(0)/(1 + zf. The evolution of the CIB 
photon number density is shown in the right panel of 
Fig. El 



APPENDIX D: SOLUTION OF THE 
BOLTZMANN EQUATIONS 



We can express the system of partial integro- 
differential equations ((6|) in terms of a system of ordinary 
integro-differential equations, 

£ i = -H£ i -b i {z,£ i ), (Dl) 
Z t = [l3 l (z,£ l )-T l (z,£ l )]Z l + (l + z)£f(z,£ i ) 1 (D2) 

where we have defined [3i{z,E) = dEbi(z,E) and 
Zi(z,E) = (1 + z)Y l (z,£ i (z,E)). The quantity £i(z,E) 
gives the energy that a particle of type i had at redshift 
z if we observe it today with energy E and take into 
account CEL. The effective source term in Eq. (|D2[) is 



Cf = d 



^2 J dE d E £j 7jj (z, £j ,Ei)- 



+ z 



(D3) 



where £j(z , E) and Zj(z, E) are subject to the boundary 
conditions £j(0,E) = E and Zj(z max ,E) = 0. The flux 
of CRs or neutrinos at z = can be expressed as 
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FIG. 6: Comparison of our results (dashed lines) for the extra-galactic proton and cosmogenic neutrino flux with previous 
calculations. Left Panel: Proton spectra from (solid lines; top to bottom) Ref. [62] (their Fig. 6 with n = 0, 7 = 2.55, 
£max = 10 12 ' 5 GeV), Ref. [| (their Fig. 14 with n = 2.4, 7 = 2.6, E max = 10 12 GeV), and Ref. [H (their Fig. 3 with 7 = 2.57, 
n = 3.3. -B max = 3 x 10 12 GeV). Right Panel: The spectra of cosmogenic v e + i/ e and + from the CMB background from 
Ref. [641 ] (their Fig. 4). We have normalized our calculations (dashed lines) to match the high energy fall-off. 



Ji(E) = —Zi(0,E) 

47T 



47T 



dz exp 



dz 



(l + z')H(z') 



H(z) 



Cf{z,£i{z,E)) 



(D4) 



In our calculation we use a logarithmic bin-size of the 
proton and neutrino energies with A log 10 E = 0.05 be- 
tween 10 5 GeV and 10 15 GeV. For the numerical evolu- 
tion of the differential equations (|D1|) and (|D2|1 we choose 
a step-size Az = 10~ 4 from to 8. The corresponding 
step-size in the propagation distance Ar = cAt is then al- 



ways smaller than the proton interaction length. In Fig. [6] 
we compare our results to previous calculations of cosmic 
proton and cosmogenic neutrino spectra, adopting the 
same values of the Hubble parameter, source evolution 
and distribution and the form of the injection spectrum. 
The agreement is satisfactory in all cases. 
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